Plateau border bulimia transition: 
discontinuities expected in three simple experiments on 2D liquid foams 
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We describe the geometry of foams squeezed between two solid plates (2D GG foams) in two main 
asymptotic regimes: fully dry floor tiles and dry pancakes. We predict an abrupt transition between 
both regimes, with a substantial change in the Plateau border radius. This should be observable 
in different types of experiments on such 2D GG foams: when foam is being progressively dried or 
wetted, when it is being squeezed further or stretched, when it coarsens through film breakage or 
Oswald ripening. 

PACS numbers: 47.20.Dr, 83.80.Iz, 47.57.Bc, 68.03.Cd 



I. INTRODUCTION 



As a physical systern, liquid foams display many inter- 
esting properties Among these, 2D foams attract 
a widespread interest : their local structure is more eas- 
ily observable than that of three-dimensional foams, and 
more amenable to theoretical calculations. 

In the present work, we present a comprehensive 
(though approximate) geometric description of 2D foams 
squeezed between two solid plates ("GG" foams). A 
thorough examination of these geometric properties sug- 
gests non-trivial consequences which should be observ- 
able through simple experimental protocols. 

This approach is extended to variational properties like 
stresses and dilatancy in a separate work [j, |5| . 

In Section [TTl we will go through an original and de- 
tailed geometrical description of a 2D GG foam. In Sec- 
tion [nil we will highlight, as a surprising consequence, 
the existence of abrupt variations of the Plateau border 
radius with tunable parameters of the foam, and describe 
three different experimental situations where such varia- 
tions should be observable. 



II. 2D GLASS-GLASS FOAM GEOMETRICAL 
VADEMECUM 

In the present section, we shall provide a geometri- 
cal description and some corresponding results for two- 
dimensional foams squeezed between two glass plates. 
More specifically, the calculations are conducted for or- 
dered^ monodisperse foams. Of course, if some local disor- 
der is present, the results should not be affected tremen- 



dously. But the foam should nevertheless be spatially 
homogeneous: there should be no large-scale gradient in 
cell thickness, bubble volume, etc. Indeed, such a gradi- 
ent would generate subtle eff'ects related to the osmotic 
pressure and force balances, which would feed the discus- 
sion substantially. 

In the remaining of this article, we call P the typi- 
cal perimeter of the bubble, defined as the perimeter of 
the rounded polygon that constitutes the bubble, as seen 
from above, and R the corresponding radius of curvature 
of the Plateau border (see Fig. [1] left) . We call H the 
distance between both solid plates, and Rps the radius of 
curvature of the pseudo Plateau borders (see Fig. [U bot- 
tom right). We also call fl the bubble volume, and finq 
the volume of liquid per bubble. In the monodisperse 
case, the liquid volume fraction (j) thus verifies: 



I.e. 



liq 



1 



(1) 
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We assume that fi remains constant: we target experi- 
mental conditions under which the applied stresses are 
not sufficient to compress the gas phase significantly. 
This is the case at atmospheric pressure, unless the bub- 
ble size is on the order of a micron or smaller. 



A. Geometry: floor tile versus pancake regime 

The main contributions to the quantity of liquid per 
bubble are pictured on Fig. [T] and can be calculated from 
simple geometrical arguments: 

l^iiq ~ (2 - ^/2) P Rl^ + -tt)R^H (2) 

The first term corresponds to the volume of the pseudo 
Plateau borders. As seen from above, they correspond 
to the light grey regions in Fig. [1] Each portion of their 
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FIG. 1: Pancake conformation of a bubble squeezed between 
two solid plates (distance H): approximate geometry. Left; 
top view. The variable P denotes the average perimeter of the 
bubbles in such a top view (outer perimeter of the light grey 
region). The variable Q is the volume of the bubble gas (full 
thickness of the white region and part of the thickness of the 
light grey region), and A is defined as ^ = Q,/H. The variable 
^tot is the total (gas and liquid) projected surface area per 
bubble (white and light grey and black regions). Right: two 
different (magnified) cross-sections (side views) with match- 
ing greyscale. As seen from above, the contact between two 
bubbles, for an ordered monodisperse foam, is typically along 
a straight line (left). In this region, the pseudo Plateau border 
(see bottom right drawing) has a uniform curvature of radius 
7?ps. By contrast, in the genuine Plateau border region (left 
drawing), the section of the gas- liquid interface is approxi- 
mately elliptical in shape (top right drawing), with principal 
radii of curvature ri and R at mid- height, and Rps and zero at 
the plates. The width Dr of the curved region is intermediate 
between Rps and H/2 while ri is larger than H/2. 



interfaces has the shape of a quarter of a circular cylinder 
(see Fig. [U bottom right). The second term in Eq. ^ 
corresponds to the genuine Plateau borders (black re- 
gions in Fig. [1]). 

Eq. ([2]) indicates that the squeezed 2D foam lies be- 
tween two asymptotic regimes depending on volume frac- 
tion and geometry. They are pictured on Fig. [21 

1. When the Plateau border radius R is much larger 
than the sample thickness H (regimes AB and CD 
of Fig. [5]), each bubble takes the form of a thick 
"pancake" , and its edge is like a half cylinder with 
radius H/2, which is the largest accessible value 
for the radius of the pseudo Plateau borders (sec 
paragraph III Bl below for a finer description). 

2. In the reverse limit, the bubbles are shaped more 
like "floor tile", with sharp edges (regime E of 
Fig. [2]): this time, the Plateau borders are like 
fine threads pinned on both solid plates, and each 
pseudo Plateau borders resembles a stretched, fine 
thread, glued on one of the solid plate and joining 
the attachment points of two Plateau borders; in 



FIG. 2: Log-log representation of the regimes of a glass-glass 
2D foam with low liquid fraction {(p <^ 1), in terms of the 
bubble perimeter P, the Plateau border radius R and the cell 
height H. 

The bubble perimeter P is measured at mid-height of the 
cell: it is the outer perimeter of the light grey ribbon in 
Fig. \T\ Such a foam can be found in two main configura- 
tions. In the fioor tile situation (regimes E and FG) the 
pseudo Plateau borders are much thinner than the cell height 
{R = Rpa <^ H). By contrast, in the pancake regime (AB 
and CD), although the overall liquid volume fraction cj> is still 
much smaller than 1, facing pseudo Plateau borders almost 
join [H — 2Rpa <^ H ^ R), and each bubble has a pancake- 
like shape. More precisely, it is useful to subdivide the pan- 
cake situation into regimes AB and CD defined by Table |T] 
in regime AB (resp. CD), the liquid content of the Plateau 
border (resp. pseudo Plateau border) dominates. The cor- 
responding expressions for the liquid volume fraction and for 
several geometrical quantities are indicated in Table HIl Note 
that regime FG, which starts with P ~ H, corresponds to 
a non-realistic 2D GG foam, where the height is larger than 
the perimeter. The limit of a truly two-dimensional foam is 
obtained with H —> oo. Note also that when P = 2n R (hori- 
zontal axis), the bubbles become independent in the liquid. 



this case, pseudo and genuine Plateau border radii 
coincide. 

On Fig. [21 wc have also pictured regime FG, which does 
not correspond to realistic 2D GG foams: it is useful to 
obtain the limit of ideal 2D foams, for which the solid 
plates are so far apart that the volume of the pseudo 
Plateau borders can be entirely neglected. This regime, 
although unrealistic, bridges the gap between 2D GG 
foam considered here and ideal 2D foams. This is useful 
in particular for discussing static dilatancy [1, H, d, 0] ■ 

In the pancake regimes (AB and CD), one has Rps — 
H/2 and hence: 

n^^^^ ~ {1/2 - tt/8) P H^ + {2V3 - tt) R^H (3) 

In regime AB, the Plateau border contribution dominates 
(black regions on Fig. [1]). Hence: 

^{2V3-tt)R^H (4) 
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Regime 



AB 



CD 



FG 



Conditions 



R<P 

{<!> < <t>.] 



W > PH 
(dominant Plateau borders) 



R-" < PH 
(dominant pseudo-Plateau borders) 



R<H 
(floor tile) 



R> H 
(pancake) 



P> H 
(2D) 



R<H 
(stretched 2D) 



P> R 

> < (Ac) 



TABLE I: Four regimes for a 2D glass-glass foam. Regimes 
AB and CD correspond to pancake-shaped bubbles, while 
regime E corresponds to a foam made of floor tile shaped bub- 
bles. In regimes CD and E, most of the liquid is located in 
the pseudo Plateau borders, while in regime AB, the Plateau 
borders themselves have a greater volume. Regime FG does 
not correspond to stable 2D GG foams (they would destabi- 
lize into 3D foams), but the limit _ff — > oo corresponds to an 
ideal, truly two-dimensional foam. 



Conversely, in regime CD, the pseudo Plateau borders 
(light grey regions on the top-view in Fig. [T]) contain 
most of the liquid: 



"liq 



(l/2-7r/8)Pi?^ 



(5) 



In regimes E and FG {R ^ H , floor tile regime), the 
radius of curvature Rps of the pseudo Plateau borders is 



equal to that of the Plateau borders: 



R 



EFG 
'ps 



R 



and Eq. ^ reduces to: 



nEFG 



[(2 - 7r/2) P + (2^/3 - tt) H] R^ 



Far away from regime FG {P ^ H ^ R), one has: 
nf,^c^i2^n/2)PR^ 



(6) 



(7) 



(8) 



The above distinction between the pancake regime and 
the floor tile regime, based solely on the overall shape of 
the bubbles (rounded edges versus sharp edges), is suffi- 
cient to provide the results of Eqs. ^ and ([7]) concerning 
the liquid volume per bubble in the various regimes (see 
Table im. 

But for other purposes, in particular concerning dila- 
tancy one must refine the geometrical description of 
the bubbles in such a GG foam in the pancake regime. 
That is the scope of the remaining of the present Section. 



B. Pancake regime with elliptical cross-section 

Only numerical tools such as Surface Evolver Q can 
provide an exact description of the bubbles when a more 
thorough geometrical description than the above approx- 
imation is needed. 



In the present paragraph devoted to the pancake 
regime, we consider an intermediate approximation based 
on an elliptical description of the pseudo Plateau border 
cross-section and derive various geometrical properties. 
A similar derivation (with a simpler geometrical descrip- 
tion) was carried out recently [9| in the context of coars- 
ening experiments. 

The shape of an inter-bubble film is a surface whose to- 
tal curvature is related to the pressure difference between 
both bubbles: 



1 

Pi 



J_ _ Pb - PA 

P2 
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(9) 



where pi and p2 are the principal radii of curvature, 
PA and pb are the gas pressure in bubbles A and B, 
and where the surface tension of the film is assumed 
to be equal to twice that of a single interface, 7. In 
the present work, for simplicity we assume that all bub- 
bles have equal pressures, so that the films are planar: 
Pi = P2 = 00. 

Since the films are assumed planar, the pseudo Plateau 
borders are circular in cross-section (radius Rps, see 
Fig. [1] bottom right) and the total curvature of the gas- 
liquid interface is therefore l/i?ps. 

In the Plateau border region, the meniscus is roughly 
toroidal in shape: the Plateau border is axisymmetric 
(radius R), and the interface is curved more strongly in 
the vertical direction (with a radius of curvature close to 
H/2). Let r(z) be the equation for its surface, in cylin- 
drical coordinates as measured from its axis of symmetry. 
The constant total curvature condition reads: 



1 



(l + r'2)3/2 y, (1 ^/2)i/g 



= const = (10) 

Rus 



where r' = dr/dz and r" — (Pr/dz^. The first term is 
the curvature within a meridian plane, while the second 
term is the curvature within a plane perpendicular to the 
meridian curve. The value of the constant is chosen so 
as to match the total curvature found for the same gas- 
liquid interface in the pseudo Plateau borders, namely 
l/i?ps. 

Of course, Eq. (jTU]) for the surface shape can be solved 
numerically. Here, we shall simply take the limit R^ H 
(which corresponds to asymptotic regimes AB and CD), 
where bubbles are not pressed very strongly against one 
another. This has two consequences: (i) the inter bubble 
films are very small in height [H — 2Rps <C H) and (ii) 
the shape of the interface in the Plateau border can be 
approximated as elliptical in cross-section. 

With this elliptical approximation (see Fig. [1] top 
right), the various geometrical parameters can be ob- 
tained very easily. At mid- height, the azimuthal radius 
of curvature is simply the Plateau border radius R. Since 
the total curvature is known to be l/i?ps, the radius of 
curvature ri in the meridian plane at mid-height can be 
obtained through: 



1 
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R 
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(11) 
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Since the length of the greater semi-axis of the eUipse is 
simply H/2, one has: 



(12) 



With these two equations and keeping in mind that R ^ 
H, one can then deduce the various lengths that appear 
on Fig. [T] (top right), including the smaller semi-axis Dr 
given by D% = Rl^ri: 
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(13) 
(14) 
(15) 



Note that the current description in terms of Plateau 
borders with an elliptical cross-section and pscudo 
Plateau borders with an circular cross-section is not en- 
tirely consistent: the crossover between both regions im- 
plies that the interface adopt an intermediate shape near 
the cuspy corners of the Plateau border. In the present 
limit R ^ H, this discrepancy can be safely neglected. 
Similarly, the elliptical approximation for the solution of 
Eq. pO|) is sufficient. Nevertheless, all subsequent results 
will not be expressed beyond the first order in H/R. 

The height e of the interbubble film is given by: 



H — 2i?ps ~ - — 



IIP 

3 "R 



(16) 



Note that with a stronger geometrical assumption (circu- 
lar rather than elliptical cross-section for the torus in the 
Plateau border), Marchalot et al. [§] obtain almost the 
same result. Only the numerical prefactor differs (1/2 in 
Ref. [i| instead of 1/3 in Eq. HID above). 



C. Liquid volume per pancake 

Let us now refine Eq. ([2]) with the elliptical approxi- 
mation presented above. 

The first term, (2 — 7r/2)Pi?pg. corresponds to the 
pseudo Plateau borders (light grey regions of Fig.[T]when 
seen from above) and already assumes that their cross- 
section is circular. Hence, this part of the equation re- 
mains correct in the refined description, provided the 
value of i?ps is taken from Eq. dH]) and provided we sub- 
tract the part of the perimeter that corresponds to the 
Plateau borders: 



H^{P - 2nR) 



2H 
~3R 



(17) 



where 27r R is the sum of the six Plateau border perime- 
ters of a typical bubble {i.e., the total length of the sep- 
aration between light grey and black regions in Fig. [1]). 



FIG. 3: Contribution to the liquid in a Plateau border. 
Within a Plateau border, the gas-liquid interface has roughly 
the shape of the outer surface of an elliptical torus (top left). 
Hence, part of the liquid in the Plateau borders of a bubble 
consists in the region between such an elliptical torus and the 
circular wall of a vertical cylinder (bottom left) with the same 
outer radius R. To calculate the volume of this contribution 
more easily, we reduce the height of both the torus (top right) 
and the cylinder (bottom right) by such a factor that the torus 
has a circular cross-section. 



The second term, (2^3 - tt) R'^H, still correctly repre- 
sents the part of the Plateau borders that spans the entire 
gap H between both plates (black regions on Fig.[T|). We 
must add, however, the volume of liquid located above 
and below the elliptical region of the interface (which 
corresponds to the curved part of the grey ribbons that 
surround the black Plateau border core regions). If we 
put together all these Plateau border edges for one single 
bubble, we obtain a fully axisymmetric region of thick- 
ness H that is described on Fig. [3] (left). Its volume 
can be deduced from that of a rounded flat cylinder with 
thickness 2 Dfi and with a toroidal edge (see Fig. [3] right). 
The volume of such a rounded flat cylinder (top right) can 
be expressed as: 



2Drtt{R-Di 



2ti{R- Dr) 



(18) 



where the flrst term is the volume of a flat cylinder with 
radius R — Dr and height 2 Dr, the second term is that 
of a semi-cylinder whose radius is Dr and whose length 
is the perimeter of the first flat cylinder. Finally, the 
third term is the correction that accounts for the cur- 
vature of the semi-cylinder, and it is equal to the vol- 
ume of a sphere of radius Dr. Subtracting these terms 
from the volume 2Dr hR^ of the outer flat cylinder (see 
Fig. [3]bottom right) and multiplying this result by a fac- 
tor H/{2Dr), wc obtain the desired liquid contribution: 



H 

20}-. 



tt{A-ti)D\R- 



7r(10-37r) 



(19) 



Combining this with Eq. (jl7p and with the second term 
of Eq. wc finally obtain: 



^^n,c.i^i/^P+(2V3-.)i?^F-i^^ 

^ 3 7r(10-37r) iJ4 
+ 12^ + 72 1^ 



(20) 
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The relative magnitude of the first three terms changes 
between regimes AB and CD (see Fig. [5]), while the fourth 
and fifth terms are always negligible in both regimes: 

regime AB : » Ti » Tg » T4 » T5 ^^i) 

regime CD : Ti > {T^ andTg) > T4 > Ts 

Because the fourth and fifth terms arc always negligible, 
we shall retain only the first three terms in Eq. ((20|) : 



^ABCD 
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1 - 



2H 



{2VS-Tr) R^H (22) 



which refines the result of Eq. ([3]). In other words, we use 
the elliptical approximation to obtain the value i?ps of the 
pseudo Plateau border radius, see Eq. pi)) . But we then 
calculate the volume as if the shape of the interface in 
the Plateau borders were identical to that in the pseudo 
Plateau borders, i.e., with circular (radius Rps) rather 
than elliptical sections (see Fig. [T] bottom right). This is 
legitimate because in the limit of small liquid fractions 
{(f) ^ 1), the total length of the Plateau borders {2tt R) 
is much smaller than that of the pscuso Plateau borders 



D. Liquid volume fraction 



From Eqs. ([7]) and (|22p . we derive the liquid volume 
fraction (/) ~ riiiq/(^tot H) hi the foam, both in the pan- 
cake regime and in the floor tile regime: 



/ABCD 



4-7r HP 



8 At. 



2H 

Jr 



+ (2^/3- 



Ato 



/EFG 



R^P 
Atot H 



(2\/3-7r) 



Ato 



(23) 
(24) 



The corresponding asymptotic values of (j) in sub- 
regimes arc indicated in Tabic [TTl It is useful to express 
the liquid fraction in terms of the volume O of the bub- 
ble itself. Using (j>/{l — (j>) = AtotH(j)/n and the total 



volume of the bubble and liquid ^tot H 
above equations become: 



n + fliin, the 



aAbcd 



4-TT H^P 



aABCD 



aefg 



aefg 



8 n 

+ {2V3-tt) 

4 - TT i?2 P 



1 - 



2H 

SR 



R^H 

+ (2V3-7r) 



R^H 

n 



(25) 
(26) 



E. Specific surface area 

With the same approximation as above, let us calculate 
the total surface area of a bubble. As it is directly related 



to the foam energy, it is useful when deriving the elastic 
modulus of the foam 0]. 

The contribution from the top and the bottom of a 
bubble to the in-plane component is the white region of 
the top view in Fig. [TJ 

2 I Aot - [2V3 -j:{R- Dnf] - {P - 2^ R) i?p,} 

In this expression, the term —2\/iR^ corresponds to one 
third of the dashed-linc triangle with edge length 2R in 
Fig. [1] and similarly for the other Plateau borders of the 
bubble. The third term, tt {R~ , corresponds to the 
dark grey sector in the triangle. The next term, — (P — 
2tt R) Rps, corresponds to the wall surface area that does 
not touch the bubble along the pseudo Plateau borders 
(light grey ribbon). 

Making the approximation Dji ~ Rps in the Plateau 
border region as mentioned at the end of Paragraph III CI 
and using Rps <^ P, this becomes: 

2Atot- {^V3~2tt)R^ -2PRps (27) 

The vertical films contribution is: 

P{H- 2Rps) (28) 

The contribution from the menisci, considered as cir- 
cular quarter cylinders with radius Rps, can be written 
as: 



2P 



Rps de^nPRps 



(29) 



Hence, the specific surface area S of the foam, equal to 
the bubble total surface area divided by the total volume 
^tot H, includes all three contributions above: 

Actios ~ 2Aot- (4\/3-27r)i?2 -2Pi?p, 
+P{H -2Rps) + Tr P Rps 



A, 



^-(4-.)^ 



Aot H 



(30) 



The final result is then derived both in the pancake 
regime where Rps is given by Eq. (fH)) and where Rps ^ R 
(regimes AB and CD of Fig. [2]) and in the floor tile regime 
where Rps = R (regimes E and FG): 
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Atot 



R^ 



Atot H 
4-7r PH 



-(4V3- 27r) 



(4-7r) 

R^ 
Aot H 



6 ^tot R 
PR 



Atot H 



(31) 



(32) 



These results are reproduced in Table HTl 
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Quantity Eqs. 


Value 


Regimes 


pseudo Plateau _Rps 
border radius ([6]ll4ll 


2 3R/ 

R 


ABCD 
EFG 


volume 
of liquid 
per bubble 


(2^3 -■k)R^ H 
^PR^ 


ABFG 
CD 
E 


liquid 
volume fraction 

(0 = :^) mm 


4-7r PH 
8 Atat 
i~TT P 
2 ^tn,, H 


ABFG 
CD 
E 


specific E 
surface area (pTll32)) 


Z V V — ^TT ) XL 

2 1 ■n--2 P 
H 2 >ltot 

H ^ Atat 

p 

Atot 


AB 
CD 

E 

FG 



TABLE II: Geometrical properties of a two-dimensional glass- 
glass foam. The numbers refer to the relevant series of equa- 
tions and the letters to the regimes of Fig. (2] pancake regime 
(AB and CD) and floor tile regime (E and FG). In each ex- 
pression, the terms are ordered by decreasing magnitudes. 



III. PREDICTED NON-TRIVIAL 
DISCONTINUITIES IN THREE 
EXPERIMENTAL CONDITIONS: THE BULIMIC 
PLATEAU BORDERS 

We shall now analyse the geometrical quantities de- 
rived in the previous section and show that they predict 
non-trivial discontinuities. 

When preparing an undeformed 2D GG foam, the ex- 
perimental conditions determine the value of three inde- 
pendent variables, for instance the liquid fraction 0, the 
bubble volume and the cell thickness H. 

Once these quantities are fixed, the bubble perimeter 
P depends principally on f2 and H, with the scaling 



n 



H, 



(33) 



since P^ scales like the bubble surface area ^tot • It also 
depends more weakly on the degree of disorder of the 
foam and on the bubble size distribution. Finally, it is 
also sensitive to the foam deformation, and this may be 
at the origin of static dilatancy, as we show elsewhere [J . 

The Plateau border radius R also results from (p, Q. 
and H and from the foam equilibrium. 

Then other quantities, such as those derived in Sec- 
tion [IT] above, are also determined. 



in regime CD, cf) does not depend on R (see Table HI]) . A 
more refined analysis, based on Eq. (PS]) . reveals that 
there is a slight dependence of on R: 



8 rj( 



A-TT mp 



which shows that 



1 



?(2^/3-7^) i?^ 2 H 



HP 3 R 



(34) 



"TT^ changes only by a factor 
close to unity when R changes from small values compa- 
rable to H to much larger values comparable to \/H P. 
In other words, the Plateau border radius R depends very 
strongly on the other parameters in this region. When 
JJ4/3 p2/3 <^R^ <^HP, Eq. (Uni) yields: 



n 



A-TT mp 



n 



Similarly, when < P < iJ^/s pi/3. 



1 _ 8 
^ 4-7r WT 



(35) 



(36) 



Thus, abrupt transitions can thus occur for the Plateau 
border radius R without any significant changes in other 
characteristic foam parameters. 

This effect can be described in physical terms as fol- 
lows. When the foam goes across regime CD, the fact 
that the Plateau border radius R is the only variable 
that varies significantly, indicates that some liquid is ex- 
changed between the Plateau borders and the pseudo 
Plateau borders. But as can be seen from Table HI in 
regime CD the volume of the Plateau borders is small 
compared to that of the pseudo Plateau borders. Hence, 
this exchange of liquid is of limited relative importance 
for the foam, which explains why other variables are af- 
fected only marginally. 



B. Three experimental situations 

Let us now imagine three experimental situations 
where the present "bulimia" effect should be apparent 
in a GG-foam. The foam should be subjected to: (1) 
progressive drying, (2) coarsening (through film rupture 
or Oswald ripening), (3) progressive squeezing. Again, 
let us insist on the fact that the present predictions rely 
on the assumption that the foam is almost at equilibrium 
and that there is no spatial inhomogeneity in the foam. 



A. Description and origin of the transition 



Drying foam 



The above calculations have a striking consequence: an 
intrinsic instability of regime CD leads to discontinuities 
in the dimension of the Plateau borders when varying 
continuously (j), H or P. 

In Section [TTl for simplicity, the liquid fraction was 
expressed in terms of P, H and or ^tot- It appears that 



Let us consider a GG-foam being progressively dried 
(decreasing 0), with constant bubble thickness H and 
volume f2, hence almost constant perimeter P. 

Starting near the liquid-solid transition, with liquid 
fraction (f> comparable to unity, the foam evolves verti- 
cally upwards in the diagram of Fig. [3^, at first in regime 
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(b) dr 




P/H 



(c) coarsening 
1 



(d) squeezing 

-1/2^2/3 1 




FIG. 4: Structural evolution of a 2D GG foam undergoing 
three different experiments. The rapid evolution through 
regime CD is symbolized by double arrows, (a) Upon drying, 
the foam goes from regime AB to regime E (in the notations of 
Fig. EJ. By contrast, when it is squeezed or when it coarsens 
spontaneously, it evolves from regime E to regime AB. (b-d) 
Schematic Plateau border radius evolution during three types 
of experiments, (b) Drying: constant P and H, decreasing (j). 
(c) Coarsening: constant cj> and H, increasing Q, and P. (d) 
Squeezing: constant (ji and $1, decreasing H and increasing P. 



AB, then in regime CD (very rapidly as a result of the 
strong variation of R mentioned in Section IIII Ap , and 
finally in regime E. 

The variations of the Plateau border radius R are plot- 
ted on Fig. [41d as a function of the liquid fraction 0. Ini- 
tially, it is comparable to the bubble perimeter, it de- 
creases as in regime AB, it then jumps down through 
regime CD to reach roughly -ff/2, then it again decreases 
as in regime E: 



Rab 



P 

pV4> 



H < RcD < VPH 
Rb ~ yfPH^ 



(37) 
(38) 
(39) 
(40) 



Of course, Eq. ([33|) must be used if the thickness H and 
the bubble volume VL are prefered to the perimeter P as 
constant parameters. 



2. Coarsening foam 

Let us now consider a GG-foam that coarsens progres- 
sively (increasing typical bubble volume O and perime- 
ter P), with constant thickness H and liquid fraction^. 
This coarsening can result from various phenomena 
When inter-bubble film are not very stable, neighbour- 
ing bubbles may coalesce and the average bubble size 



increases. Alternatively, coarsening can result from con- 
tinuous Oswald ripening: gas diffusion between neigh- 
bouring bubbles leads to a net flux from small bubbles 
to large bubbles as a result of the larger pressure in the 
smaller bubbles; large bubbles grow even larger, small 
bubbles vanish and the size distribution evolves to larger 
length scales. 

Starting near the 2D stability transition, with a hori- 
zontal size (or perimeter P) comparable to the thickness 
H , the foam evolves as shown on the diagram of Fig. |4^, 
at first in regime E, then in regime CD (very rapidly as 
a result of the strong variation of R mentioned in Sec- 
tion |TTT^, and finally horizontally in regime AB. 

The variations of the Plateau border radius R are plot- 
ted on Fig. 01: as a function of the bubble volume Vl. It 
increases like fJ^/* in regime E, it then jumps up through 
regime CD from H/2 to roughly H/ y/Tfi^ then it again 
increases like -v/fi in regime AB: 



i?E ^ i^lHf'^^ 
H < RcD<H/^ 
Rab ^ <t>/H 



(41) 
(42) 
(43) 



Of course, Eq. ([55)) must be used if the perimeter P is 
prefered to the bubble volume Q. as the main variable. 

In coarsening experiments such as those conducted by 
Marchalot et al. a (delicate) measurement of the 
Plateau border radius R will be necessary to observe the 
phenomenon described above clearly. 



3. Squeezing a foam 

Let us now consider a GG-foam that is being squeezed 
progressively (decreasing thickness H) at constant bub- 
ble volume CI and liquid fraction (j). 

Starting near the 2D stability transition, with a hori- 
zontal size (or perimeter P) comparable to the thickness 
H, the foam evolves as shown on the diagram of Fig. 2^, 
at first in regime E, then in regime CD (very rapidly as 
a result of the strong variation of R mentioned in Sec- 
tion |TTT^, and finally horizontally in regime AB. 

The variations of the Plateau border radius R are plot- 
ted on Fig. [Hi as a function of the sample thickness H. 
It is non-monotonic: it decreases like H^^^ in regime E, 
it then jumps up through regime CD, then it increases 
like 1/VH in regime AB: 

i?E ~ (nH)^^^^ (44) 

f7l/3 ^2/3 < < ^1/3 ^1/6 (45) 

Rab ^ V^HH (46) 

Of course, Eq. ([55[) must be used if the perimeter P is 
prefered to the thickness H as the main variable. 
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IV. CONCLUSION 

By examining the geometry of 2D "GG" -foams (lo- 
cated between two parallel solid plates), we found that 
these foams should display a rather strong "bulimia" 
transition where Plateau borders swallow (going from re- 
time E to regime AB) a large amount of liquid when the 
foam parameters arc varied in a narrow region. This 
effect, which is specific to 2D "GG"-foams, should be 
observable in at least three types of experiments, for 
which we provided first-order predictions. Of course, 
there may be other experimental situations. Let us in- 
sist again on the fact that the present predictions are 
restricted to foams almost at equilibrium and with no 



spatial inhomogencity. This "bulimia" effect should af- 
fect other phenomena, such as dilatancy; this is discussed 
elsewhere [3, 0] ■ 
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